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ABSTRACT 

fN) ' We present new high-resolution spectroscopic observations (A/AA « 40000) of 18 

bright northern Cepheids carried out at David Dunlap Observatory, in 1997. The 
measurements mainly extend those of presented in Paper I (Kiss 1998) adding three 
more stars (AW Per, SV Vul, T Mon). The spectra were obtained in the yellow-red 
spectral region in the interval of 5900 A and 6660 A, including strong lines of sodium D 
and Ha. New radial velocities determined with the cross-correlation technique and the 
bisector technique are presented. The new data are compared with those recently pub- 
lished by several groups. We found systematic differences between the spectroscopic 
Q\ \ and CORAVEL-type measurements as large as 1-3 km s _1 in certain phases. 

Q\ . We performed Baade-Wesselink analysis for CK Cam discovered by the Hipparcos 

satellite. The resulting radius is 31±1 R Q , which is in very good agreement with recent 
period-radius relation by Gieren et al. (1999). It is shown that the systematic velocity 
differences do not affect the Baade-Wesselink radius more than 1% for CK Cam. 

Observational pieces of evidence of possible velocity gradient affecting the individ- 
ual line profiles are studied. The full-width at half minimum (FWHM) of the metallic 
lines, similarly to the velocity differences, shows a very characteristic phase depen- 
dence, illustrating the effect of global compression in the atmosphere. The smallest 
line widths always occur around the maximal radius, while the largest FWHM is as- 
sociated with the velocity reversal before the minimal radius. Three first overtone 
pulsators do not follow the general trend: the largest FWHM in SU Cas and SZ Tau 
■ occurs after the smallest radius, during the expansion, while in V1334 Cyg there are 

only barely visible FWHM- variations. The possibility of a bright yellow companion 
of V1334 Cyg is briefly discussed. The observed line profile asymmetries exceed the 
values predicted with a simple projection effect by a factor of 2-3. This could be asso- 
ciated with the velocity gradient, which is also supported by the differences between 
individual line velocities of different excitation potentials. 
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1 INTRODUCTION 

Detailed understanding of Cepheid variables, as strictly pe- 
riodic pulsating stars with supersonic atmospheric motions, 
requires spectroscopic observations in a wide wavelength in- 

* Based on observations obtained at David Dunlap Observatory, 
Canada 
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terval with high spectral resolution and good phase cover- 
age. The published high-resolution (at least A/AA ~ 20000— 
30000) Cepheid spectroscopy addressed to the kinematic and 
dynamic phenomena covers only a sample of few bright stars. 
The recently presented extensive optical observations are as 
follows: Wallerstein et al. (1992) - W Sgr, k Pav, S Mus, 
S Nor, (3 Dor, Y Oph, U Car; Butler (1993) - FF Aql, 
8 Cep, r) Aql, X Cyg; Breitfellner & Gillet (1993a,b,c) - 
S Cep, j] Aql, S Sge, X Cyg; Sabbey et al. (1995) - 5 Cep, 
T) Aql, C Gem, X Sgr; Baldry et al. (1997) - I Car; Butler 
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& Bell (1997) - S Cep, rj Aql, X Cyg; Gillet et al. (1999) - 
S Cep. Obviously, the overwhelming majority of the north- 
ern classical Cepheids - except a few (3-4) stars - suffers 
from the observational neglect in this point of view. This 
fact turned our attention to the topic and the first results 
based on echelle spectroscopy were already published in Kiss 
(1998). 

The problem of spectral line formation in a moving stel- 
lar atmosphere was studied theoretically by Karp (1978). He 
showed that the different Doppler-shifts observed in lines of 
different strengths depend on the ionization balance varia- 
tions as well as on the velocity gradient. Butler (1993) pre- 
sented phase dependent velocity differences between lines 
of different excitation and ionizations. He found up to 5 
km s~ velocity residuals compared to high excitation po- 
tential (EP) Fe I lines during the phase of rapidly decreasing 
velocity. This important issue was explored in much details 
by Sabbey et al. (1995). They discussed the importance of 
line profile asymmetries when measuring radial velocities. 
In an earlier paper of Wallerstein et al. (1992) a comparison 
between different measuring techniques was done with simi- 
lar purposes. One of the main conclusions in these papers is 
that assigning a Doppler-shift to a specific portion of the line 
profiles is a difficult and ambiguous task. Recently, Vinko et 
al. (1999) studied the limitations of precise radial velocity 
measurements in Cepheid atmospheres and concluded that 
there is a natural limit of a few km s _1 , which could not 
be overstepped. They also showed that CORAVEL-type ra- 
dial velocity (i.e. cross-correlation velocity smoothed over 
the whole visible spectral region) curves may differ from 
those obtained from a selected set of spectral lines. 

Another question concerning the variable line profiles 
is the atmospheric motions in classical Cepheids described 
by the turbulence variations. Turbulence is considered to 
be one of the most important line broadening mechanisms 
besides rotation, although its physical meaning is still un- 
clear. CORAVEL-type measurements (Benz & Mayor 1982, 
Bersier & Burki 1996) led to the determination of turbulence 
variations in certain Cepheids, although these studies have 
neglected the effect of possible velocity gradient. Bersier & 
Burki (1996) pointed out the similarity of acceleration and 
turbulence curves with local maxima around the velocity 
reversal. Breitfellner & Gillet (1993a, b,c) used a nonlinear, 
nonadiabatic model of pulsation in order to determine an 
excess FWHM of the observed spectral lines. They analyzed 
four stars with different periods and found different turbu- 
lent velocity curves. Recently, Gillet et al. (1999) presented 
FWHM variations for S Cep with excellent phase coverage, 
concluding that the strongest peak in the turbulent velocity 
curve at <^>=0.85 is associated with the global atmospheric 
compression. They also suggested the possibility of turbu- 
lence amplification induced by shock waves of moderate in- 
tensities. 

The main aim of this paper is to present new spectro- 
scopic observations for a larger sample of Cepheid variables 
covering a wider period range. The (above mentioned) ear- 
lier papers dealt almost exclusively with fundamental pul- 
sators, while in the short-period domain there are a few 
bright overtone Cepheids (e.g. SU Cas, DT Cyg, SZ Tau). 
Bersier & Burki (1996) suggested the different turbulence 
behaviour of classical and s-Cepheids as a possible mode- 
discriminator. Therefore, beside the new radial velocity 



Table 1. The list of programme stars. Periods and epochs were 
improved in five stars (typesetted in bold face) using photometric 
data published in Paper I. The epoch and period for T Mon were 
taken from Evans et al. (1999). The remaining ephemcridcs arc 
from Szabados (1991) and Paper I (unchanged). 
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data, we examine the fine spectral differences between fun- 
damental and overtone pulsators. 

The paper is organised as follows. The observations are 
described in Sect. 2. New radial velocities together with the 
earlier ones are discussed in Sect. 3, Sect. 4 deals with the 
line profile analysis (variations of the FWHM and the asym- 
metries along the pulsational cycle). A summary is given in 
Sect. 5. 



2 OBSERVATIONS 

The high-resolution spectroscopic observations were carried 
out at David Dunlap Observatory with the echelle spectro- 
graph attached to the 74-inch telescope in 1997 Septem- 
ber. The detector and the setup was the same as used by 
Kiss (1998), except the cross-grating, which was replaced by 
the one with 600 lines/mm. Therefore, we detected 12 or- 
ders with dispersion of 3.8 A mm -1 giving a resolution of 
A/AA « 40000 at Ha. All orders covered about 40 A be- 
tween 5860 and 6660 A with ^25 A gaps between them. In 
this spectral region we could identify more than 100 lines, 
the strong sodium D doublet, Ha and photospheric metal 
(mostly Fe I) lines. One order was highly, another four were 
partially contaminated by atmospheric telluric lines, which 
were monitored every night with observing the rapidly ro- 
tating bright telluric standard HD 177724 (V=2.99 mag, 
spectral type A0V). As telluric lines change with airmass 
and time, we did not apply telluric line corrections, only ex- 
cluded the affected spectral region from the further analysis. 
The exposure times were between 5 and 40 min depending 
on the target brightness and atmospheric conditions in order 
to reach a signal-to-noise ratio of 50-250. The considerably 
high cosmic ray contamination after 30 minutes did not al- 
low the use of longer exposures. 
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RT Aur, 1997.09.26/27. 
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Figure 1. A sample spectrum for RT Aurigac. 



The list of observed stars is presented in Table 1 . It con- 
tains mostly the same stars as in Paper I. SV Vul, AW Per 
and T Mon was observed only in 1997, therefore, they were 
not present in the spectroscopic programme of Paper I. We 
continued the spectroscopic monitoring of two double-mode 
Cepheids, CO Aur and TU Cas. 

The spectra were reduced with standard IRAF tasks 



including bias removal, flat-fielding, cosmic ray elimination, 
echelle orders extraction (with the task doecslit) and wave- 
length calibration. For this calibration, we obtained two 
ThAr spectral lamp exposures, immediately before and af- 
ter every stellar exposure. We applied careful linear interpo- 
lation between the spectral lamp exposures, as slow wave- 
length shifts in the obtained spectra were observed caused 
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Figure 2. Cross-correlation velocities from two spectral orders 
for HD 187691 during the two observing seasons. The systematic 
shift due to an instrumental problem (see text) is 0.50 km s -1 . 



by the movement of telescope. Each order was normalized 
to the continuum by fitting cubic splines. The regions of 
strong lines (sodium D, Ha) were omitted from the contin- 
uum fitting. A typical spectrum with 12 orders for RT Aur 
is shown in Fig. 1, while the list of the marked spectral 
lines is presented in Table 2. Note, that spectral lines in 
long-period Cepheids (i.e. stars with lower effective temper- 
atures) are more numerous and stronger, thus Table 2 gives 
only a rough overview of the obvious spectral features. 



3 RADIAL VELOCITIES 

Radial velocities presented in Table 3 were determined with 
two different methods. In order to compare the data with 
those of in Paper I, we obtained a set of radial velocities 
by cross-correlating two spectral regions of the Cepheid and 
IAU standard star spectra with the IRAF task fxcor. We 
chose HD 187691 (spectral type F8V, v rad = +0.1 ± 0.3 
km s _1 ) as an overall template star. The cross-correlated 
regions were those of 6188-6220 and 6405-6435 A, where 19 
photospheric lines were identified (see Table 2). 

The stability of radial velocity measurements was tested 
by the standard star itself. We plotted the cross-correlation 
velocities in Fig. 2, where a systematic shift of 0.50 km s _1 
is present between data obtained in 1996 and 1997. The 
necessary shift for matching the data is the same for the dif- 
ferent wavelength regions. Furthermore, other non-variable 
stars (a Per, HD 22484) showed similar shift. This turned 
out to be an instrumental effect, probably the entrance slit 
of the spectrograph was not uniformly illuminated by the 
spectral lamp. Therefore, a systematic asymmetry occured 
in the emission line profiles of the ThAr lamp introducing a 
systematic subpixel-shift in the wavelength calibration. We 
determined the mean correction of 0.5 km s _1 and applied 
it to the calculated radial velocities. In order to guarantee 
the homogeneity of the data reduction procedure, we have 
also re-reduced our spectra obtained in 1996. The newly 
calculated radial velocities do not differ from values in Kiss 
(1998) more than 0.1-0.14 km s _1 , the mean deviation is 
about —0.05 km s _1 . This is less by a factor of 3 than the 



Table 2. List of identified strong (I CO re < 0.9) spectral lines. 
Lines used for bisector velocity determination are typesetted in 
bold face. Rest-wavelengths and excitation potentials were taken 
from Moore, Minnaert & Houtgast (1966). 
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Table 3. The observed heliocentric radial velocities. V\: cross-correlation velocities; V2: mean bisector velocities. The spectra of double- 
mode Ccpheids TU Cas and CO Aur did not allow determination of accurate single-line velocities due to the increased observational 
scatter caused by the faintness of the targets. 
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Table 3. (cont.) The observed heliocentric radial velocities. 
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internal accuracy of measurements estimated to be about 
±0.15 km s" 1 . All data are tabulated in Table 3. 

Three stars (FF Aql, V1334 Cyg and S Sge) show re- 
markable changes of the 7-velocity due to orbital motion in 
binary sytems. The best example with good phase coverage 
is V1334 Cyg, where V 7 changed about 11 km s _1 in one 
year. This is almost a half of the full orbital velocity am- 
plitude, which is 27.8 km s _1 , according to Evans (1995). 
The sytematic velocity changes in FF Aql and S Sge are 
2.3 km s _1 and 27.7 km s _1 , respectively. These stars have 
well-determined orbits (see Evans 1995), thus we did not try 
to re-calculate their orbital elements. 



3.1 Comparison of velocity curves 

It is interesting to compare our new cross-correlation veloc- 
ity curves with those obtained by others with other instru- 
ments for the same stars, especially with the CORAVEL- 
type velocities (Bersier et al. 1994), because those have been 
extensively used in many studies of Cepheid variables. This 
comparison can be seen in Fig. 3 for T Vul, 5 Cep and X Cyg 
where the continuous lines mean CORAVEL-velocities and 
symbols denote others' data including ours. These three 
stars have well-covered velocity curves and they are used 
in order to represent shorter and longer period Cepheids. In 
the case of 5 Cep data of Breitfellner & Gillet (1993a) and 
Butler (1993) were also plotted. The velocities of Breitfellner 
& Gillet are based on a few selected Fe I lines, while But- 
ler's data have been obtained with the very precise iodine- 
cell technique. Note, that Butler (1993) reported that his 
data have no absolute zero-point, and a 2-3 km s _1 shift is 
necessary to match his S Cep-velocities with those of others. 
We also added a 3 km s _1 correction to these velocities and 
found perfect agreement with our ccf-velocity curve as well 
as the velocities of Breitfellner & Gillet (1993a). Contrary 
to these, there is a significant difference between all these 



velocities and the CORAVEL data on the ascending branch 
of the velocity curve of 8 Cep. 

An earlier version of this graph appeared in Paper 
I showing larger differences between the datasets. Those 
larger shifts were partly due to a slightly longer period used 
in Paper I for 5 Cep. In this paper we used an updated, 
shorter period (see Table 1) that eliminated the larger dis- 
persion of the phased velocities on the ascending branch. 
However, a shift between the CORAVEL and other types of 
velocities is still present. The same can be observed in the 
case of X Cyg and T Vul, although for T Vul the largest 
differences occur during the velocity reversal. As in the pre- 
vious case, phase shifts due to slightly incorrect epochs or 
periods have been eliminated by carefully taking into ac- 
count the period variations of these stars, therefore all the 
differences can be considered real. Moreover, Vinko et al. 
(1998) found similar deviations between digital ccf- and 
CORAVEL-velocities of Type II Cepheids. 

All of these comparisons are consistent with each other 
and suggest that for Cepheid variables the CORAVEL-type 
velocities systematically differ from those obtained by direct 
spectral line measurements. Possible reasons for the cause 
of this discrepancy are the well-known effects that perturb 
the line profiles of pulsating variables, namely the velocity 
differences and/or the line asymmetries. These effects are 
phase-dependent, so they may be capable of explaining why 
the deviations of the CORAVEL-velocities get stronger at 
certain phases. On the other hand, it is a bit surprising that 
our digital cross correlation (ccf-) velocities (and also those 
of Vinko et al., 1998) agree much better with line bisector- 
or iodine cell-velocities than with CORAVEL-type ones. The 
digital cross-correlation velocities should be (and probably 
they actually are) also sensitive to the line profile distur- 
bances mentioned above. It is possible that another system- 
atic effects are also present in the CORAVEL-type measure- 
ments (perhaps associated with the alignment of the mask 
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3.2 Line level effects and velocity differences 
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Figure 3. Comparison of recently published radial velocities for 
T Vul, <5 Cep and X Cyg. 



or the usage of much more spectral features in a much wider 
wavelength interval) that contribute to this disagreement. 
The velocity data obtained with CORAVEL have very good 
phase coverage and also very good inner precision, therefore 
they contain very relevant information on pulsating stars. 
Therefore we conclude that the CORAVEL data should be 
compared with other velocity curves obtained with other 
techniques before using them for deriving astrophysical in- 
formation that may be sensitive to a few km s _1 systematic 
differences. 



Single-line radial velocities were also calculated for a set 
of well-defined, unblended lines observed in both seasons. 
We used the line bisector technique (see Wallerstein et al. 
1992 for the definition) to derive radial velocities measur- 
ing Doppler-shifts of unblended Fe I, Ni I, Ti 1 and Si II 
lines in the 6100-6450 A region. We chose the 0.7 bisector, 
because our tests showed that it gives the most accurate ve- 
locities. The precision of the individual line velocities were 
determined using artificial digital spectra with random noise 
co-added. The precision turned out to be ±500m s _1 and 
±150m s- 1 for spectra of S/N=50 and S/N=200, respec- 
tively (Vinko et al. 1999). Mean radial velocities were cal- 
culated and plotted in Figs. 4-5 with the corresponding ccf 
data (Vi) for eight stars covering a wide period range (1.9-27 
days). Radial velocity differences between lines of different 
excitation potentials (EP) and the mean metallic velocities 
(in sense of Vii ne minus V mC an) are also plotted in the sub- 
panels of Figs. 4-5. The referee suggested to also include Ha 
velocities, however, the Ha profile is strongly affected by var- 
ious atmospheric phenomena and, consequently, it is quite 
misleading to associate one velocity to the distorted line pro- 
file (see, e.g. Wallerstein et al. 1992, Vinko et al. 1998). A 
detailed analysis of Ha-observations, especially concerning 
the Ha emission in long-period Cepheids, will be presented 
in a forthcoming paper (Vinko & Kiss, in preparation) . Here 
we want to focus on the velocity behaviour of metallic lines. 

A few important conclusions can be drawn based on 
Figs. 4-5. First, the average 0.7 bisector velocities are in 
good agreement with the ccf data at level of 0.5-1 km s _1 . 
This is an expected result, since the cross-correlation tech- 
nique results in a mean Doppler-shift for lines in the selected 
wavelength region. Second, there are clear examples of line 
level effect, especially for low-EP (< 2 eV) and high-EP (> 6 
eV) lines. Our results are very similar to those of presented 
by Butler (1993) for 4 stars. The Si II 6347.1 A line is a 
very good indicator of the strength of level effects. Its devi- 
ations, however, cannot be explained by kinematic velocity 
gradients (Butler et al. 1996), because the 7- velocity of the 
Si Il-velocities is quite different from the 7-velocity of the 
average velocity curve. 

It is visible in the lower panels of Figs. 4-5 that the 7- 
velocities of some other lines also differs from the 7- velocity 
of the average radial velocity curve. This was also reported 
by Butler (1993), and it gives further support to the con- 
clusion of Sabbey et al. (1995) that path conservation may 
not be valid for the integral of the observed velocity curves 
of Cepheids over the whole pulsational cycle. 

The level effects get stronger with increasing period 
starting from 1-2 km s _1 for P < 7 days up to 6-8 km s _1 
for P > 7 days. Note that the ±2 km s _1 scattering in the 
velocity differences of SU Cas and T Vul at later phases is 
due to observational uncertainties rather than real physical 
effects. In the spectra of these stars the lowest and highest 
EP lines become very weak at these phases which degrades 
the precision of the velocity measurement. In the case of the 
other stars the observed velocity differences can be consid- 
ered significant, especially for the longer period ones. Un- 
fortunately, the data of T Mon have very bad phase cover- 
age, but the observed tendencies suggest similar behaviour 
to X Cyg, namely the low-EP lines have positive velocity 
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Figure 4. Radial velocity data for SU Cas, T Vul, <5 Ccp and r\ Aql. Upper panels show the cross-correlation (open circles) and mean 
bisector velocities (pluses). Differences between the individual and mean bisector data for selected lines (pluses - 6358.6879 A, Fc I 0.86 
cV; open squares - 6327.604 A, Ni I 1.68 cV; solid down triangles - 6411.658 A, Fc I 3.65 cV; crosses - 6408.026 A, Fe I 3.69 eV; solid 
circles - 6347.095 A, Si II 8.12 cV) are plotted in the lower panels. A radial velocity offset of 2 km s was applied for clarity. 



residuals while the high-EP lines have large negative veloc- 
ity residuals in the 0.6-1.0 phase interval. 

We tried to give an order of magnitude estimate of the 
dependence of the velocity gradients on the stellar param- 
eters from the observed velocity differences. This requires 
the knowledge of the geometric distance between the line 
forming regions of low- and high-EP lines, which is diffi- 
cult to determine in a dynamic atmosphere. As a very crude 
first approximation, we assumed that the lowest and high- 



est E.P. lines form at the top and at the bottom of the 
photosphere, and that the width of the photosphere is in- 
versely proportional to the gravity, thus Ax p ~ ctg 1 where 
q is roughly constant (Gray, 1992). The absolute value of 
the velocity gradient is then V« w Av/Ax v w pAv r g/ct 
where p is the conversion factor of radial to pulsational ve- 
locities. We chose Av r = t> r (FeI 0.86 eV) -iv(SiII 8.12 eV) 
and a constant p-factor of p =1.36. Therefore we obtain 
log Vi> = logAiv + logp + constant. It is visible that the 
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Figure 5. Radial velocity data for S Sge, £ Gem, X Cyg and T Mon. The symbols used are the same as in Fig. 4. 



velocity gradient has weaker dependence on the pulsational 
period or the radius of the Cepheid, because while Av r in- 
creases with increasing period (or radius via the period- 
radius relation), the gravity decreases at the same time. 
Thus, the velocity gradients are expected to be roughly the 
same for short- and long-period Cepheids, provided the sim- 
ple scaling with gravity is indeed valid. Using the gravity 
values presented in Paper II (Kiss & Szatmary 1998) it can 
be seen that the velocity gradients calculated in this way 
have the same order of magnitude in the case of S Cep and 
X Cyg. The absolute value can be estimated assuming fur- 
ther that the thickness of the Cepheid photospheres can be 



approximated with the scaling of the width of the solar pho- 
tosphere (about 500 km) to log g — 1 — 2, resulting in a value 
of Vw = 10~ 5 s _1 . However, in the case of strong shock 
waves forming in the photosphere, the velocity gradient can 
be much larger than this, at least temporarily. 

3.3 Baade-Wesselink radius of CK Cam 

The Baade-Wesselink analysis is a widely used tool for ra- 
dius determination of pulsating variables. Here we do not 
want to discuss the recent theoretical and observational ef- 
forts which have been done to clarify the problems of this 
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Figure 6. Baade-Wcsselink analysis of CK Cam. 

method, since it is beyond the scope of present paper. We re- 
call e.g. Gautschy (1987) for a general review. On the other 
hand, recent surface-brightness techniques (e.g. Gieren et al. 
1997) are based on different infrared colours, which enables 
a reliable separation of temperature and radius variation. 

This subsection is addressed to CK Cam, which was 
discovered by the Hipparcos satellite (Makarov et al. 1994). 
Because of the relatively new discovery of this star, there 
are only few available photometric observations. Thus, the 
sophisticated Baade-Wesselink implemantations or surface- 
brightness methods cannot be used. That is why we per- 
formed a simpler analysis following the original assumptions 
of Wesselink (1946). This approach was used recently by Ba- 
log et al. (1997) for Type II Cepheids, and it is based on BV 
photometry and a full radial velocity curve. 

We have taken the published BV photometry of Berd- 
nikov et al. (1996). In order to reduce the numerical uncer- 
tainties, we have fitted low-order (3-5) Fourier-polynomials 
to the light, colour and radial velocity curves. They are plot- 
ted in Fig. 6, where the displacement curve is also shown. 
The inferred radius is 31±1 R©. This is in very good agree- 
ment with the prediction of recent period-radius relation of 
Gieren et al. (1999) giving for P=3.29 days a radius of 31.5 
R,0. This relatively large radius implies that CK Cam is a 
regular Type I Cepheid (i.e. young supergiant star), which 
is also supported by the low galactic latitude (8.7°). 

The effect of phase-dependent velocity differences was 
found to be much lower than the estimated systematic errors 
due to the initial simplifications. By replacing Vi to V2 data 
the resulting radius changed only 0.2 Rq (0.6%). 



4 LINE PROFILE VARIATIONS 

Different line broadening mechanisms occuring in a pulsat- 
ing atmosphere can be globally described by the variation of 
the Full Width at Half Minimum (FWHM). Recently, Gillet 
et al. (1999) have used this parameter to trace the turbu- 
lent velocity variations in 8 Cep. From nonlinear, nonadia- 
batic pulsational models they conclude that the main factor 
governing the line broadening processes is the global com- 
pression/expansion of the atmosphere. Shock wave effects 
turned out to be much weaker in their models. 

In order to get an overall picture about the period de- 



pendence of this phenomenon, we have determined FWHM 
curves for 14 stars similar to that of Gillet et al. (1999). 
We chose several unblended lines of different excitation po- 
tentials (EP) and found the shape of FWHM curves not 
depending on the actual EP value. We plotted the resulting 
curves for 14 stars in Fig. 7 calculated for line A 6411.658 
A, which is a neutral iron line of EP=3.65 eV. 

The general trends are very similar for the majority of 
the stars. The smallest FWHM always occurs very close to 
the phase of maximal radius. On the other hand, the largest 
FWHM is usually associated with the velocity reversal point 
(around <f> w 0.8-0.85). This is just before the smallest ra- 
dius, when the global compression is the strongest, domi- 
nating the line broadening. 

There are three deviating stars in our sample. Two of 
them (SU Cas and SZ Tau) have the smallest line width 
in similar phases than in other stars, however, the largest 
FWHM occurs later, between 0=0.95-1.00 (see Fig. 8). 
These variables are thought to be first overtone pulsators 
(referred as s-Cepheids) which could be a likely explana- 
tion for the deviation. Bersier & Burki (1996) suggested a 
separation between classical and s-Cepheids based on their 
radial velocity and turbulence variations. This aspect of at- 
mospheric phenomena, i.e. the effects of overtone pulsation, 
has been rarely studied theoretically. Evans et al. (1998) 
explained the strong period variation of the first overtone 
Cepheid Polaris by the complexities of the envelope acoustic 
cavity. Nevertheless, there is a number of pieces of evidence 
for a possible connection between the atmospheric turbu- 
lence phenomena and the mode of pulsation. In our sample, 
there are two other s-Cepheids (DT Cyg and FF Aql) , which 
do not follow the suggested distinction, therefore, we cannot 
draw a firm conclusion on this issue. 

The third deviating star is V1334 Cyg, which is a well- 
known binary Cepheid with a bright companion. This star 
neither follows the overall pattern of FWHM- variations (i.e. 
largest FWHM occurs around 0=0.8-0.9), nor fits the ex- 
pectations on the absolute value of the full- width. As can be 
seen in Fig. 7, its line profile is quite stable around 0.67 A, 
while there is a slight decrease around 4> ~0.75. A possible 
explanation for these peculiar behaviour could be the effect 
of a bright companion which contributes significantly to the 
observed line profiles (see below). 

The asymmetry of spectral lines of pulsating stars have 
been extensively investigated by many authors, recently Sas- 
selov & Lester (1990), Sabbey et al. (1995) and Albrow & 
Cottrell (1994). The asymmetry parameter (AP) was defined 
in several forms in these papers. In the followings we adopt 
the form given by Sasselov & Lester (1990), namely the log- 
arithm of left-to-right half-width ratio of a given spectral 
line (in this case the Fe I A 6411.658 A). We checked the de- 
pendence of AP on choosing the spectral line and found no 
significant variation in the spectral region covered by our ob- 
servations. The precision of the AP values was tested using 
some telluric lines which are expected to be close to AP=0. 
The uncertainty of AP was found to be 1-2 percent. 

The behaviour of AP during the pulsational cycle can be 
observed in the right side of Fig. 7 for the programme stars. 
In the case of SU Cas, DT Cyg and SZ Tau only very weak 
variations can be seen, while there are much larger changes 
in the case of long-period Cepheids. The natural explana- 
tion is the smoother pulsation in first overtone stars, where 
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Figure 7. FWHM and line profile asymmetry variations of the iron line A 6411.658 A in 14 Cephcids. 



the amplitudes of pulsation are smaller and the dynamic 
effects are weaker. VI 334 Cyg again deviates from other 
low-amplitude pulsators showing asymmetry variation that 
is comparable to long-period Cepheids. On the other hand, 
the large asymmetries occuring in long-period stars are as- 
sociated with large FWHM variations, while there is no such 
phenomenon in V1334 Cyg. 

Generally, the AP variation curves resemble the radial 
velocity curves indicating that the projection effect (the cen- 
ter of the stellar disk shows maximal Doppler-shift, while the 
limb has no Doppler-shift) contributes at least partly to the 
observed asymmetries. On the other hand, the model calcu- 
lations by Sabbey et al. (1995) showed that there are other 
factors that also cause line asymmetries (mainly the varying 
depth of line formation over the pulsation cycle) producing 
larger observed asymmetries. It is visible in Fig. 7, that the 
amount of asymmetry variation (10-20 percent) is about the 
same for all fundamental pulsators in our sample. These new 
observations provide additional confirmation of the theoret- 
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gest a possible effect of the mode of pulsation on the atmospheric 
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Figure 9. Line doubling in V1334 Cyg compared with similar 
Cepheids in similar pulsational phases. 



ical results obtained by Sabbey et al. (1995) and Albrow & 
Cottrell (1994). 



4.1 The companion of V1334 Cygni 

The low-amplitude pulsator V1334 Cyg is unique among 
other low-amplitude Cepheids in many respects, as was 
shown in the previous sections. Here we try to construct 
a self-consistent explanation of all these peculiarities by as- 
suming a bright yellow companion star moving in a close 
orbit around V1334 Cyg A. Both the unusually constant 
FWHM and large AP variation can be explained with the 
presence of the spectral lines of the component star with 
slightly Doppler-shifted to the lines of V1334 Cyg A. The 
secondary lines make the combined line profiles broader, 
thus, the variation of FWHM due to pulsation is obscured. 
Furthermore, the Doppler-shift of the pulsation displace the 
Cepheid lines from the companion lines producing higher 
asymmetries. In parallel with these effects, the measured 
amplitude of the radial velocity variation decreases because 



the secondary lines move the centers of the combined lines 
closer to their equilibrium positions. All these effects can be 
seen in a numerical model using simulated data that is a 
subject of a follow-up publication (Kiss, in prep.). 

We made a few very simple considerations on the pos- 
sible nature of the hypothetic component star. The visual 
amplitude of VI 334 Cyg is about 0.2 mag being too small 
compared to other short-period Cepheids. If we consider 
V1334 Cyg to be a fundamental mode pulsator, its am- 
plitude should be about 0.4-0.6 mag (see, e.g., CK Cam, 
RT Aur, T Vul). Assuming a virtual amplitude decrease 
due to the significant amount of secondary light, one can 
compute the luminosity ratio L(2)/L(Cepheid). In order 
to get an amplitude of 0.2 mag instead of 0.5, a ratio of 
L(2)/L(Cepheid)~3 is required suggesting an unphysical 
yellow supergiant, which is too bright and still undetected. 
However, a more consistent result is given by the assumption 
of first overtone pulsation. In that case the intrinsic visual 
amplitude is about 0.4 mag (see, e.g., SU Cas, DT Cyg, 
SZ Tau, FF Aql), and the observed 0.2 mag implies a ratio 
L(2)/L(Cepheid)«l. This is still a very bright component, 
but it gives more consistent picture of the system. If it has 
a similar spectral type than V1334 Cyg A does, than it is 
almost invisible in the IUE spectra and consequently, its de- 
tection is much more difficult in the ultraviolet. The bright 
blue companion (B7.0V) cannot be responsible for all these 
effects because it is much fainter than the Cepheid itself in 
this spectral region and also it is a fast rotator with v sin i 
around 200 km s _1 (Evans, personal communication). 

We found some observations directly suggesting the 
presence of a bright secondary component in the line profiles. 
Fig. 9. shows line profiles of V1334 Cyg at certain phases 
corresponding to the highest asymmetries (thick lines). As 
a comparison, line profiles of other Cepheids with similar pe- 
riod at the same phases are also plotted (thin lines). It can 
be seen that a) the lines of V1334 Cyg are much broader 
and b) a multiple structure is resolved in V1334 Cyg that is 
absent in the other two Cepheids (note, that FF Aql is also 
a binary but its component is much fainter to be detected in 
this spectral range). Both the broad lines and the multiple 
structure strongly support the hypothesis of a bright yellow 
companion. 

We did not observe systematic shift of the FWHM and 
AP curves between 1996 and 1997, while during this period 
the 7-velocity of V1334 Cyg changed by 11 km s _1 due to 
orbital motion. Consequently, the component that causes 
the orbital motion cannot be the same that is responsible 
for the line profile anomalies. We propose that there is a 
close, yellow component of V1334 Cyg A on an orbit with 
low inclination producing no observable Doppler-shifts, but 
affecting the line profiles of the pulsating star. These two 
stars move together along an orbit around a massive third 
star with an amplitude of 27.8 km s _1 (Evans 1995). Because 
the blue component detected by the IUE may not be that 
massive third star, the system of V1334 Cyg might be a 
complicated multiple system. 



5 SUMMARY 

The results presented in this paper can be summarized as 
follows. 
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1. We made high-resolution echelle spectroscopy for 18 
northern Cepheids in the yellow-red spectral region, between 
5900 A and 6660 A. New radial velocity data with internal 
accuracy of about 0.15 km s _1 were calculated using the 
digital cross-correlation technique. Single-line radial veloci- 
ties were also determined by the bisector technique and their 
average values were compared with the ccf-data. There is a 
good agreement at a level of 0.5-1 km s _1 between these 
two techniques applied for the same observations. 

2. We compared our new measurements with recently 
published spectroscopic and CORAVEL-type velocities. We 
found i) perfect agreement with conventional spectroscopic 
data and it) systematic differences as large as 1-3 km s _1 
between CORAVEL and our data in certain phases. Pos- 
sible explanations for this discrepancy are the line profile 
perturbing effects of velocity gradient and line asymmetries, 
although our digital ccf-data are less affected than those of 
obtained by the CORAVEL technique. 

3. Line level effects were studied for a set of spectral 
lines with different excitation potential. We found very clear 
examples of line level effect for low-EP (<2 eV) and high- 
EP (>6 cV) lines. The Si II 6347.095 A line is a very good 
indicator of the strength of level effects. The observed level 
effect get stronger with the increasing period, starting from 
1-2 km s _1 for P<7 days up to 6-8 km s _1 for P>7 days. 
Using these velocity differences we estimated the velocity 
gradients and their dependence on the stellar parameters. 
Our approximative approach gave a roughly constant value 
of Vv = KT 5 s _1 . 

4. We performed a Baade-Wesselink analysis for 
CK Cam discovered by the Hipparcos satellite. The inferred 
radius is 31±1 R suggesting a regular Type I Cepheid. 

5. The phase dependent behaviour of FWMH and asym- 
metry variations were examined. The fundamental and over- 
tone pulsators seem to follow different trends. In the case 
of fundamental pulsation, the smallest FWHM always oc- 
curs very close to the phase of maximal radius, while the 
largest FWHM is associated with the global compression of 
the atmosphere. In the first overtone Cepheids the largest 
FWHM is shifted toward the temperature maximum, around 
(f> = 0.95 — 1.00. The observed line asymmetries are con- 
sistent with theoretical predictions involving the effects of 
projection and varying depth of line formation over the pul- 
sation cycle. 

6. We present observational pieces of evidence for a 
bright, yellow companion of V1334 Cyg based on the fine 
spectral peculiarities of this Cepheid. This suspected de- 
creases the measurable light and radial velocity amplitudes 
by a factor of 2 and may explain the other peculiarities of 
the line profile variations. 
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